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NUMERICAL SOLUTION OF THE HYPERSONIC 


VISCOUS-SHOCK-LAYER EQUATIONS FOR LAMINAR, TRANSITIONAL, 

AND TURBULENT FLOWS OF A PERFECT GAS OVER 
BLUNT AXIALLY SYMMETRIC BODIES 

E. C. Anderson* and James N. Moss 
Langley Research Center 

SUMMARY 

The viscous-shock-layer equations for hypersonic laminar, transitional, and turbu- 
lent flows of a perfect gas over two-dimensional plane or axially symmetric blunt bodies 
are developed and presented. The resulting equations are transformed to a shock-body 
oriented coordinate system and solved by an implicit finite-difference scheme. The eddy 
viscosity is approximated by a two-layer model which has been used extensively in 
turbulent- boundary -layer analyses. Methods for defining the boundary-layer-thickness 
parameters (required in the eddy-viscosity model) which are consistent with boundary- 
layer theory are presented. 

Numerical solutions of the equations of motion for turbulent viscous shock layers 
are compared with those of a turbulent-boundary-layer analysis which is not corrected 
for the entropy variation of the inviscid flow outside the boundary layer. The two solu- 
tion procedures are in good agreement in the region of flow where the effects of stream- 
line swallowing by the boundary layer are negligible. For the downstream regions, where 
streamline swallowing effects are present, the expected differences in the two solution 
procedures are evident. However, to establish the accuracy of the downstream solutions, 
comparisons with a boundary -layer analysis which accounts for the effects of streamline 
swallowing are necessary. 

INTRODUCTION 

The nose blunting required to reduce the surface heat-transfer rates to an accept- 
able level in the stagnation region of a hypersonic vehicle results in strong entropy gra- 
dients through the boundary layer in the downstream region. In the stagnation region, 
the vorticity interaction between the boundary layer and the external inviscid flow must 
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be considered in the analysis if the Reynolds number is sufficiently low. Mayne and 
Adams (ref. 1) present detailed comparisons of a number of boundary-layer analyses with 
the viscous-shock-layer analysis developed by Davis (ref. 2). The results presented by 
Mayne and Adams show that the viscous -shock-layer analysis correctly accounts for 
viscous-inviscid vorticity interaction, streamline swallowing, and viscous -induced pres- 
sure effects. Although the viscous-shock-layer analysis was found to be the more accu- 
rate solution technique, Mayne and Adams demonstrate that a conventional boundary-layer 
analysis corrected for streamline swallowing is sufficiently accurate and is computation- 
ally more efficient than the viscous-shock-layer analysis. The comparisons presented 
by Mayne and Adams are restricted to laminar flows of a perfect gas. 

As a result of the increased boundary -layer thickness in turbulent flows, the effects 
of viscous -induced pressure and the strong entropy gradients produced by the bow shock 
are more pronounced. For turbulent flows, the solution of the viscous-shock-layer equa- 
tions should provide a substantially more accurate result than would be obtained with a 
boundary -layer analysis corrected for streamline swallowing without additional correc- 
tions to include interaction of the boundary layer with the external inviscid flow. 

For flow conditions in which heat transfer from the shock layer by radiation must 
be considered, the formulation of the viscous-shock-layer equations provides an impor- 
tant simplification. The presence of the radiative heat flux term in both the inviscid and 
viscous governing equations necessitates an iterative solution procedure of both sets of 
equations to achieve the required coupling. Since the necessary coupling of the inviscid- 
and viscous-flow equations is contained within the viscous-shock-layer equations, the 
solution of this set of governing equations should be more accurate and computationally 
more efficient than an iterative solution technique which uses a boundary-layer-type 
analysis. However, applications of the viscous-shock-layer equations to turbulent flows 
have not been reported. An independent investigation of the viscous -shock-layer equa- 
tions for turbulent flow has been developed by Eaton and Larson (ref. 3) concurrently with 
the present study. The results presented by Eaton and Larson are restricted to the thin- 
shock-layer assumption and to pointed bodies, as a result of their definition for the loca- 
tion of the boundary -layer edge. The present study considers the complete viscous-shock- 
layer equations for turbulent flow and presents definitions for the boundary -layer thickness 
and incompressible boundary-layer displacement thickness which are consistent with 
boundary -layer theory. Both the present study and that of reference 3 are applicable to a 
perfect gas. Either analysis can, with minor modifications, be applied to equilibrium 
flows. 

It is noted that the boundary-layer-thickness parameters are introduced into the 
solution by the two-layer eddy-viscosity model used in the present study and would be 
unnecessary if an invariant-turbulence model were available. The eddy viscosity is 
approximated by use of the formulation proposed by Cebeci (ref. 4). 
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Numerical solutions of the viscous-shock-layer equations for turbulent flow are 
compared with solutions from the boundary -layer analysis developed by Anderson and 
Lewis (ref. 5). The two solution procedures show acceptable agreement over the portion 
of the body where the effects of streamline swallowing by the boundary layer are negligi- 
ble, and the expected differences are obtained in the region far downstream of the nose. 
The boundary -layer analysis does not account for streamline swallowing. 

SYMBOLS 


A + 


damping factor (eqs. (23) and (24)) 


a 0 

a l 


quantity defined by equation (A9e) 


quantity defined by equation (AlOf) 


a 2 

C 


f 


quantity defined by equation (AlOg) 


skin-friction coefficient 


H 


specific heat at constant pressure 

u 2 

defined quantity, h + — 


H t total enthalpy, H + 


h specific enthalpy, hy 


j flow index: 0 for plane flow; 1 for axisymmetric flow 


k thermal conductivity 


lCp eddy thermal conductivity 


St mixing length (eq. (21)) 


M Mach number 


* * / * 

Np r Prandtl number, p c pA 
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N Pr,T 

N Re 

n 

n + 

P + 

P 


r 

r 

s 

s 
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* 

n 



u 

U T 

V 

v+ 

X 


if j if 

turbulent Prandtl number, c pPqy ^ 

Reynolds number, P«,U*rJy£i* 
coordinate measured normal to body, n*)/r* 
normal coordinate (eq. (22)) 
pressure-gradient parameter (eq. (25)) 
pressure, p 

wall heat-transfer rate 

radius measured from axis of symmetry to point on body surface, r 
nose radius 

coordinate measured along body surface, s*/!-* 

s at beginning of transition 

transition damping factor (eq. (35)) 

temperature, TyT* e f 

temperature, 

free-stream velocity 

velocity component tangent to body surface, u^/u^ 
friction velocity (eq. (27)) 

velocity component normal to body surface, v^U^ 
scaled mean velocity component (eq. (26)), v^u^. 
transition length scale 
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a 

“ 1 ^ 2 ^ 

a 3> a 4j 


shock angle defined in figure 1 
coefficients in equation (A8) 


0 

y 


6 

5 k 

e + 

* 

V 

e 


K 


P 


M'j' 




P 


o 


angle defined in figure 1 

ratio of specific heats 

normal intermittency factor (eq. (30)) 

streamwise transition intermittency factor (eq. (34)) 

boundary -layer thickness 

incompressible displacement thickness (eq. (29)) 

normalized eddy viscosity (eq. (4)), P T j\i 

eddy viscosity, inner law (eq. (20)) 

eddy viscosity, outer law (eq. (28)) 

transformed n-coordinate, n/n s 

body angle defined in figure 1 

body curvature 

molecular viscosity, P/P*( T ref) 
eddy viscosity (eq. (1)) 

coordinate measured along body surface, | = s 
density of mixture, P*/p^o 


Reynolds number parameter, 


M*(Tref) 

n* TJ* r* 

^oo u co A n 


1/2 



<p quantity defined by equation (11a) 

X vorticity Reynolds number (eq. (33)) 

X c vorticity Reynolds number at beginning of transition 

if/ quantity defined by equation (AlOe) 

Superscripts: 

j 0 for plane flow; 1 for axisymmetric flow 

_ time-averaged quantity 

* dimensional quantity 

' total differential or fluctuating component 

" shock-oriented velocity component (see fig. 1) 

Subscripts: 

e boundary -layer edge 

s shock 

w wall 

°o free stream 

Abbreviations: 

BL boundary layer 

VSL viscous shock layer 
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ANALYSIS 


Governing Equations 

The equations of motion for turbulent viscous shock layers are derived by methods 
analogous to those presented by Dorrance (ref. 6) for the turbulent-boundary-layer equa- 
tions. The resulting equations for the shock-body coordinate system shown in figure 1 
are presented in nondimensional form. For turbulent flow, time-averaged quantities are 
implied by the nomenclature, and the laminar-flow equations are obtained by neglecting 
the turbulent eddy viscosity and eddy thermal conductivity. 

The eddy viscosity is expressed in time-averaged fluctuating components as 


P t - 


(pv)V 

2 3u 
a * 

3n 


( 1 ) 


and the eddy thermal conductivity is 


k T = - 


(pv) h' 


3h 

3n 


( 2 ) 


From the definition of the static turbulent Prandtl number 

_ Vt 
Npr > T " IT 


(3) 


and the definition 

p-j. 

M 


c+ - 


(4) 


the governing equations for the flow in laminar or turbulent viscous shock layers can be 
expressed as 


Continuity: 


•^[(r + n cos 0) j pu] + JL[(1 + nx)(r + n cos 0) jpvj = 0 


(5) 
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s-momentum: 


_u + v — + UV/C 

1 + nK 3s 3n 1 + nK 


1 + Tin 3s 


= o 2 / — 


3n 


3n 


MU K 
1 + nK 


2k j cos e 
+ 


1 + nK r + n cos 6, 


" (1 + <+) r-, i 

3n 1 + nK 


MUK 


( 6 ) 


n-momentum: 


p _JL_8v + v lx . 
\1 + nK 3s 3n 


U 2 K ^ 
1 + n k) 


3n 


(?) 


Energy: 


P_JL_M + vM 

\l + nK 3s 3n / 


3p Pu z vk o J 8 

v — + = 

3n 1 + nK ] 3n 


M 


N 


N 


1 + e + 


Pr \ 3H 


Pr 


N Pr jX / 3n 


+ c p 


N, 


±™LL- Jf_i +e+ J^M + ^ 

yl + nK r + n cos OJ Np r y N Pr,T/ 9n 


State: 

p = ^^pT 

y 

In these equations, the molecular viscosity as given by Sutherland's law is 

= lj_C t 3/2 
T + C 


( 8 ) 

( 9 ) 


(10) 


where 


<p = 


P 

Npr 



1 + 


e + N 


Pr 


N 


Pr,T 


( N Pr,- 
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1 + nK 


dla) 


8 



(lib) 


c - 


(y - iJm^t* 


and 


C* = 110.33 K (11c) 

The governing equations have a hyperbolic-parabolic nature (ref. 2), the hyperbolic 
nature arising from the normal-momentum equation. If the thin- shock-layer approxima- 
tion is made, the normal-momentum equation becomes 

£P = P^JL (12) 

9n 1 + n« 

When equation (7) is replaced by equation (12), the resulting set of equations is parabolic. 
Consequently, the equations can be solved by using numerical methods similar to those 
used in solving boundary -layer problems. Equation (12) is used for an initial iteration; 
then for the final flow-field solution, equation (12) is replaced by equation (7), so that the 
thin-shock-layer approximation is removed. 

Boundary Conditions 

Conditions at the body surface .- The no-slip boundary conditions are used in this 
study. The surface conditions for tj - 0 are 

u = v = 0 (13) 

and for this study the temperature and enthalpy at the wall are specified as 

T w = H w = Constant (14) 

Conditions at the shock .- The conditions imposed at the shock are calculated by 
using the Rankine-Hugoniot relations. The nondimensional shock relations are as follows 

Mass: 


PgVg = -sin a (15) 

u” = cos a (16) 

b 
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Momentum: 


p s = + sin 2 a jl 

rM M V 



( 17 ) 


Energy: 


T = 
s 


2yM 2 sin^ a - (y - l)||(y - 1)M 2 sin 2 a + 2 


(y - l)(y + l) 2 M 0O 4 sin 2 a 


(18) 


State: 


yp q 

p = s — 

s (y - 1)T S 


(19) 


A transformation is applied to the previous nondimensional viscous-shock-layer 
equations and boundary conditions to simplify the numerical computations. The transfor- 
mation relations and the transformed equations and boundary conditions are given in the 
appendix. 


Eddy-Viscosity Approximations 

A two-layer eddy -viscosity model consisting of an inner law based upon Prandtl's 
mixing-length concept and the Clauser-Klebanoff expression (based on refs. 7 and 8) for 
the outer law is used in the present investigation. This model, introduced by Cebeci 
(ref. 4), assumes that the inner law is applicable for the flow from the wall outward to the 
location where the eddy viscosity given by the inner law is equal to that of the outer law. 
The outer law is then assumed applicable for the remainder of the viscous layer. It is 
noted that the eddy viscosity degenerates to approximately zero in the inviscid portion of 
the shock layer. The degeneracy is expressed in terms of the normal intermittency fac- 
tor given by Klebanoff (ref. 7). The expressions used in the present investigation are 
given in the following sections. 

Inner-eddy-viscosity approximation .- Prandtl's mixing-length concept is stated in 
nondimensional variables as 


.+ _ Pi 2 9u 
* a 2 ^ 9n 


( 20 ) 
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The mixing length i is evaluated by using Van Driest's proposal (ref. 9) stated as 


a 


kin! 



where 


n 


+ 


np 

^W| 

^u\ 

a 2 p 

P 

\dnj. 


1/2 


( 21 ) 


( 22 ) 


Here, kj is the Von Karman constant, which is assumed to have a value of 0.4, and A + 
is a damping factor. 

Cebeci (ref. 4) suggests that for flows with a pressure gradient, the damping factor 
be expressed as 


A + = 26(1 - 11.8P+)* 1 / 2 


(23) 


and for flows with both a pressure gradient and mass injection, 

r + r 1/2 

A + = 26 / - P-[exp(U.8v+) - l] + exp(11.8v+) \ 


where 


(24) 


and 


P + = -a 2 M _J£ 


8s 'e p 2 u T 3 


(25) 




(26) 


u T - q 


^w/9u 
P \3n / 


'w 


1/2 


(27) 


Outer-eddy-viscosity approximation . - For the outer region of the viscous layer the 
eddy viscosity is approximated by the Clauser-Klebanoff expression 
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(28) 


+ 



k 2 pu eVi,n 


cr 2 u 


where 
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k 



k 2 = 0.0168 


(29) 


and 




(30) 


Equation (30) is Cebeci's approximation (ref. 4) of the error-function definition presented 
by Klebanoff (ref. 7). 

The boundary -layer thickness 6 is assumed to be the value of n at the point 

where 


— !- = 0.995 (31) 

\ H t- 

and is determined by linear interpolation in an array of local total enthalpies. This 
definition is approximately equivalent to the usual boundary -layer definition 

-S- - 0.995 (32) 

u e 

/ 

where u e is the local value for the undisturbed inviscid flow outside the boundary layer. 

The values of the parameters kj and k 2 in equations (21) and (28) depend on the 
flow conditions being considered, as does the constant represented by the value 26 in 
equations (23) and (24). The values given are used for convenience in developing the 
numerical method. For a discussion of these parameters, see reference 10. 

Transition Models 

Provisions for both continuous and instantaneous transition from laminar to turbu- 
lent flow have been included in the present investigation. Instantaneous transition is 
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initiated when the local Reynolds number or momentum-thickness Reynolds number 
exceeds a preselected value. Continuous transition is effected by defining a streamwise 
transition intermittency factor y. . which modifies the composite eddy viscosity e + 

over an interval x. The factor y. . is initially set to zero and is evaluated when the 

IjS 

vorticity Reynolds number (proposed by Rouse, ref. 11, as the stability index) 


_ n^p 9u 
cr^p 9n 

exceeds a critical value y c . For the evaluation of y^ |, the relation 
y. y - 1 - exp(-0.412s) 

1 J S 


(33) 


(34) 


was used, where 


- 4 ( s - s o) 

s = - 

s o(* - !) 


(35) 


This model was developed by Dhawan and Narasimha (ref. 12) on the basis of the experi- 
mental data presented by Owen (ref. 13). Approximate values of x c and x are 


2000 = x c = 4000 


and 


x =2 

The values of x c and x are strongly influenced by the body shape and flow conditions 
being considered and are more appropriately defined by comparison with experimental 
data. A discussion of these parameters is given by Harris (ref. 10). 

Method of Solution 

The procedure for solving the viscous-shock-layer equations is presented herein. 
First, the finite-difference expressions used to transform the differential equations to 
algebraic equations are presented. Then the solution procedure is discussed. 

Finite-difference expressions .- The derivatives are converted to finite -difference 
form by using Taylor’s series expansions. A variable grid spacing (fig. 2) is used in the 
p-direction so that the grid spacing can be made small in the region of large gradients. 
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Three-point differences are used in the 77 -direction, and two-point, fully implicit differ- 
ences are used in the ([-direction. Truncation terms of order A£ m (first-order accu- 
racy) and either Ari n Ari n _ j or At? - &ri n _ 7 (second-order accuracy) are neglected. A 

typical finite-difference expansion of the standard differential equation (see eq. (A8)) 
gives 


■\i^m,n-l + ^n^m,n + ^n^ m ,n+l ~ ^n (36) 

The coefficients A, B, C, and D are used to represent the coefficients after the 
finite-difference expansion of equation (A8). The subscript n denotes the grid points 
along a line normal to the body surface, whereas the subscript m denotes the grid sta- 
tions along the body surface. Equation (36) along with the boundary conditions constitutes 
a system of the tridiagonal form, for which efficient computational procedures are avail- 
able. (See ref. 14.) 

Overall solution procedure .- For specified free-stream conditions and body geom- 
etry, a stagnation-streamline solution is obtained. With the stagnation- streamline solu- 
tion providing the initial conditions, the conditions at the shock providing the outer bound- 
ary conditions, and the conditions at the wall taken as the inner boundary conditions, the 
numerical solution is marched downstream to the desired body location £. The first 
solution pass provides only an approximate flow-field solution, because the following 
assumptions are used in the first solution pass: 

(a) The thin-shock-layer form of the n-momentum equation (eq. (A12b)) is used. 

(b) The stagnation-streamline solution is independent of downstream influence 
(approximation of local similarity where n 2 s = 0). 

(c) The term dn s /d([ is equated to zero at each body station. 

(d) The shock angle a is assumed to be the same as the body angle 9. 

These assumptions are then removed by making one or more additional solution 
passes. For the current study, a total of two solution passes are used since the two 
passes resulted in a converged flow-field solution. For the second solution pass, the 
thin- shock-layer form of the normal-momentum equation (eq. (A12b)) is replaced by equa- 
tion (A12a). The v component of velocity that is used in equation (A12a) is the value 
from the previous solution pass. Also, once the first solution pass has been computed, 
the values of ^ s and dn g /di[ are calculated and used in the second solution pass to 
remove approximations (b), (c), and (d). Hence, the viscous-shock-layer equations are 
solved as parabolic equations, and yet retain effects which are elliptic and hyperbolic in 
nature. This solution procedure is programed for the Control Data Corporation 6600 com- 
puter. The execution time is approximately 0.03 second per grid point for a converged 
solution. (This includes all local iterations and solution passes.) 
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Shock solution .- The shock solution procedure at any location is identical for the 
first and subsequent solution passes. However, the shock angle a is defined differently 
for the first and subsequent solution passes and is set equal to the local body angle 6 for 
the first solution pass. For subsequent solution passes, the shock angle is defined as 

a = 8 + tan 1 (37) 

1 + /Olg 

Solution procedure at station m . - The viscous -shock-layer equations are solved at 
any body station m (see fig. 2) in the order shown in figure 3. The governing equations 
are uncoupled and the dependent variables are solved one at a time, also in the order 
shown in figure 3. First, the shock conditions are calculated to establish the outer bound- 
ary conditions. Then the converged profiles at station m - 1 are used as the initial 
guess for the profiles at station m. The solution is then iterated locally until conver- 
gence is achieved. For the stagnation streamline (m = 1), guess values for the profiles 
are used to start the solution. 

The first-order equations are numerically integrated by means of the trapezoid 
rule. Each of the second-order partial differential equations is individually integrated 
numerically by using the tridiagonal formalism (eq. (36)). The global -continuity equa- 
tion is used to obtain both the shock standoff distance and the v components of velocity. 
By integrating equation (All) between the limits of rj = 0 and rj = 1 at station m, an 
implicit equation for n s is obtained. For the v component of velocity at rj, equa- 
tion (All) is integrated with respect to q between the limits of 0 to 77 . The pressure 
is determined at station m by integrating the normal- momentum equations (A12) with 
respect to rj between the limits of 1 to rj. The equation of state is used to determine 
the density. 

DISCUSSION AND RESULTS 

Comparisons are made between the present viscous-shock-layer analysis and 
turbulent-boundary-layer solutions obtained with the computer program described in ref- 
erence 15. For the data presented, fully developed turbulent flow without mass injection 
has been assumed. However, provisions have been included for the analysis of continuous 
or instantaneous transition and for mass injection at the surface. 

Results are presented for a Mach 19 flow about a 45° (total angle) hyperboloid with 
a nose radius of 0.3048 meter for free-stream Reynolds numbers of 3.1 x 10^ and 
12.4 x 10 6 . The wall temperature is taken to be one-tenth of the stagnation temperature, 
and the free-stream temperature is 256 K. The Prandtl number and the turbulent Prandtl 
number are assumed to be 0.72 and 0.90, respectively. 
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The combination of high Reynolds numbers and low-wall-temperature results in 
boundary-layer thicknesses which are less than 10 percent of the total shock-layer thick- 
ness and imposes severe resolution requirements upon the numerical method. The solu- 
tions presented have been obtained with 150 grid points across the shock layer, the ratio 
of adjacent step sizes being a constant value equal to 1.08 (geometric progression). This 
choice of the grid-point distribution results in approximately 80 points within the pre- 
dominantly viscous portion of the shock layer. More recent calculations have shown that 
100 grid points across the shock layer having the same constant step-size ratio provide 
adequate resolution. The differences in the solutions obtained with 150 and 100 grid 
points are less than ±1 percent. However, the use of 75 grid points across the shock 
layer results in solutions which differ by as much as 30 to 40 percent for the wall meas- 
urable quantities. For laminar flows with low Reynolds number and a higher ratio of 
wall temperature to stagnation temperature, as few as 50 grid points can be used. 

Distributions of boundary -layer thickness, incompressible displacement thickness, 
heat transfer, and skin friction calculated by the present analysis and by boundary -layer 
theory are shown in figures 4, 5, 6, and 7, respectively. Inspection of these results shows 
that maximum differences between the two methods are less than ±10 percent for all 
parameters. The heat-transfer distribution corresponding to a Newtonian surface pres- 
sure distribution as input to the boundary -layer computer program is shown in figure 6 
for a free- stream Reynolds number of 3.1 x 10®. Less severe conditions than those for 
the cases presented result in essentially identical solutions from the two methods of 
analysis. 

Representative comparisons of the eddy-viscosity profiles for a free-stream 
Reynolds number of 3.1 x 10® are shown for body stations Sy^r* of 0.5 and 3.5 in fig- 
ure 8. At the downstream location, maximum differences in the eddy viscosities are less 
than ±10 percent. For the body station in the nose region, differences as large as 30 per- 
cent are noted for n*yr* near the boundary-layer edge. These differences were found 
to have little influence upon the solution. Transition criteria indicate that the flow at the 
station in the nose region is laminar. 

Figures 9 and 10 show velocity profiles calculated by the viscous-shock-layer and 
boundary -layer solution methods for Nj^ e = 3.1 x 10® at s*/r* = 2 and s*y/r* = 80, 
respectively. The data shown in figure 9 for s*/r* = 2 correspond to a location where 
the effects of streamline swallowing are negligible. For this body station, the profiles 
corresponding to boundary -layer and viscous-shock-layer solutions differ by less than 
±5 percent for n*/r* = 0.86. At the downstream body station s*/r* = 80 (fig. 10), the 
agreement between the profiles calculated by the two methods is poor, since the boundary- 
layer analysis does not account for the effect of streamline swallowing. The edge condi- 
tions for the boundary-layer solution correspond to an isentropic expansion of the gas 
from the stagnation point of the body. For the viscous-shock-layer solution, the edge 
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conditions correspond to the flow which crosses the weaker portions of the bow shock, 
where the entropy jump across the shock is much less than at the stagnation streamline. 
For the 45° (total angle) hyperboloid considered, the asymptotic limit approaches the 
solution for an equivalent cone. For the body station s/'r* = 80, the conditions at the 
edge of the boundary layer as calculated with the viscous- shock-layer analysis differ 
from the cone solution by less than 1 percent. The surface pressure corresponding to 
the viscous-shock-layer solution at this station is approximately 15 percent higher than 
the corresponding inviscid cone solution. However, as previously mentioned, in the outer 
inviscid region the solution is approximately the same as the asymptotic cone solution. 

The increase in pressure through the predominantly viscous layer was also noted by 
Mayne and Adams (ref. 1) for laminar flow. For laminar flows, the viscous- shock-layer 
solution used in the present analysis is the same as that used in reference 1. Results for 
laminar flows are not presented. 

Heat-transfer distributions corresponding to the viscous -shock- layer and boundary- 
layer solutions for Np^ e = 3.1 x 10® are shown in figure 11. The two solutions show the 
expected differences, and the viscous-shock-layer analysis appears to account correctly 
for streamline swallowing. However, the boundary-layer analysis used in the present 
study does not account for streamline swallowing, and the accuracy of the viscous-shock- 
layer analysis must be assessed by comparison with an appropriately modified boundary- 
layer analysis. 

For the solution of the turbulent-viscous-shock-layer equations, the use of 
Klebanoff's intermittency factor (ref. 7) is essential, whereas in boundary-layer solutions, 
the intermittency factor may be assumed unity without significantly influencing the results 
obtained. (See ref. 16.) For the viscous-shock-layer equations, the use of a unit inter- 
mittency factor results in an increase in heat transfer and skin friction of 30 to 50 percent 
for syr* > 2. This difference in behavior of the viscous-shock-layer and boundary- 
layer solutions is the result of the nonvanishing tangential velocity gradients in the 
normal-coordinate direction. (See fig. 9.) For the boundary-layer equations, the bound- 
ary conditions imposed ensure that the gradients in the normal- coordinate direction 
approach zero at the boundary -layer edge. 

CONCLUSIONS 

Equations describing the turbulent viscous shock layers over blunt axially symmet- 
ric bodies of analytic shape are presented for hypersonic flow of a perfect gas. A two- 
layer eddy-viscosity model consisting of an inner law based upon Prandtl’s mixing-length 
concept and the Clauser-Klebanoff expression for the outer law is used in the present 
study. Methods for defining the boundary -layer thickness and incompressible boundary- 
layer displacement thickness which are consistent with boundary -layer theory are pre- 
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sented. An implicit finite -difference technique for solving the equations is given. Com- 
parisons of the present results with previously reported boundary -layer solutions are 
made for a Mach 19 flow over a 45° hyperboloid. The free-stream Reynolds numbers 
are 3.1 x 10® and 12.4 x 10®, whereas the wall temperature is one-tenth of the stagnation 
temperature. 

Results of the study lead to the following conclusions: 

1. The present results are in good agreement with classical boundary-layer results 
in regions where the effects of streamline swallowing are negligible. However, for the 
downstream locations, where streamline swallowing effects are present, expected differ- 
ences between the present results and the classical boundary -layer results are evident. 

2. For the solution of the turbulent-viscous-shock-layer equations, the use of 
Klebanoff's normal intermittency factor is essential. This is in marked contrast with 
turbulent-boundary-layer solutions, where the normal intermittency factor may be 
assumed unity without significantly influencing the results. 

3. For the flow conditions considered, a boundary-layer thickness (required in the 
eddy-viscosity model) based on a total -enthalpy ratio provides good agreement with 
boundary -layer results. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 21, 1975. 
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APPENDIX 


TRANSFORMED VISCOUS- SHOCK- LAYER EQUATIONS 

This appendix presents the transformed viscous-shock-layer equations and boundary 
conditions. First, the relations defining the transformed variables and coordinates are 
given. Next, the general equations and boundary conditions are given. Then, the special 
form of the equations for the stagnation streamline are developed along with the stagna- 
tion shock relations. 


Transformation Relations 

To simplify the numerical computations, a transformation is applied to the viscous- 
shock-layer equations. This transformation is accomplished by normalizing most of the 
variables by their local shock values. When the normal coordinate is normalized with 
respect to the local shock standoff distance, a constant number of finite -difference grid 
points between the body and shock are used. This procedure eliminates the need for 
interpolation to determine shock shape and the addition of grid points in the normal direc- 
tion as the computation moves downstream. 

The transformed variables are 


(Al) 


where bars over the quantities denote normalization with respect to the shock value. The 
transformations relating the differential quantities are 



d_ = &_ % _ 8 _ 

8s 0£ n s 77 8 n 


(A2) 


19 



APPENDIX - Continued 


where 

dn c 
n' = — - 
s d£ 

_ 9 _ = 1_B_ 
3n n g 9r? 

and 


(A3) 


(A4) 


9 2 _ 1 9 2 

c 2 2 Q 2 

3n n g 9?? 

The transformations used to express the shock-oriented velocities 
terms of the body-oriented coordinate system (fig. 1) are 


(A5) 


u" and v" in 
s s 


and 


Uc = u" sin (a + 0) + v" cos (a + 0) 

® S b 


(A6) 


v s = -Ug' cos (a + 1 3) + Vg sin ( a + 0) 


(A7) 


Transformed Equations 

After the governing equations are written in transformed variables and coordinates, 
the second-order partial differential equations are written in the following form: 


9 2 W 
9t? 2 


+ a 1 — + Q!<jW + ao + a, — = 0 
1 3rj 2 3 4 9£ 


(A8) 


The quantity W represents u in the s- momentum equation and H in the energy equa- 
tion. The coefficients to a ^ are written as follows: 


s-momentum, W = u: 


a . 


1 8a 0 


n s K 


a n 9?? 1 + n s T]K 


2 - 


a 0,s a 0/ 


jn g cos 9 


n s^s u s n s 


PUT? 


n sP s v s pv 


r + n s r? cos 6 a 2 a 0>s (l + n s 7 ?/c) a 0 a2a 0,s a 0 


(A9a) 


20 



APPENDIX — Continued 


a 2 = ~ 


Kn s 1 8a 0 n s * | 

K j COS 9 \ 

| ^s ** 

1 + n g T7K a.Q 9tj 1 + n s 7]K' 

yl + n s r )K r + n s 7] cos 8, 

^ a 0,s a 0 


n 0 2 u' 
s s s 


up 


n s Vs K 


pv 


o^a^gU + tV?*) c^a^Cl + n s ? 7 /c) a 0 


a o = - ■ 


Ps n s 2 Zap PsP_^8^ 


ff^O^iigii+ngJiKlao^ p s n s 


d V 


« 4 



2 

u s n s pu - 
(1 + n s T]K) a 0 


where 


a Q = p(l + e + ) 
Energy, W = H: 


1 9a l 

a , = + n c 


1 a.! 


j< + j cos e + 


p s n s | 

?n' s u s riup 

g2a l ,s a l 

yl + n s ? JK 


- V s pv 


a 2 = 


a 4 H s 

H c 


“3 = 


CtA = ~ 


"S 2 


a l,s a l H s 

n s 977 

"s 2 

p s u s^ 


i cos 6 , , 

+ 1 1// 


n s 8rj 1 1 + n s rjK r + n s r/ cos 8 


°\\s^ + n s^) a l 


(A9b) 


(A9c) 


(A9d) 


(A9e) 


(A 10a) 


(A 10b) 


(AlOc) 


(AlOd) 
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APPENDIX — Continued 


where 


a 2 ~ a 

^ = u “ul-ii 


2~r2,s a 2 aa 


n s 8 77 1 + T7n s K 


S MU 


(AlOe) 


. M /, . £+N Pr 


a, = — 1 — 1 + — 
1 Nn. N 


Pr V Pr,T/ 


(AlOf) 


2 N 


M 


Pr 


e + N 


N n - 1 + — 
Pr N 


Pr 


Pr,T 


N Pr ,T " 1 


. (AlOg) 


The preceding energy equation is for the thin-shock-layer approximation. When equa- 
tion (7) is used for the n-momentum equation, the following term must be added to 
equation (A 10c): 


n s v s v 


, 35 n s u s p s K -21 

[p s -ti - pu^ 1 


° 2 a ljS H s a 1 \ S 1 + 7 l n s K 


(AlOh) 


The remaining equations are written as follows: 


Global continuity: 


a 

H 


n s (r + n s r7 cos 0 ) J P u s pu 


+ — ( (r + n S T 7 cos 0) ] 


(1 + n s T/K)p v s jov - n' 77P u s pu 


= 0 


(All) 


n-momentum: 



1 + n s 77K\ v s 


V 8v\ + p^ 0V _ ^s K - -2 + p s 3p _ Q 
n s 3 t] I u s n s 0TJ • v s 1 + n s 7]K P s u s n s v s 0T J 


(A12a) 
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APPENDIX - Continued 


which, if the thin-shock-layer approximation is made, becomes 

8p Ws 2 * --2 
dri p g (l + n s r?K) 


(A12b) 


State: 


p = pT (A13) 

Equations (A8) to (A13) along with the appropriate boundary conditions are the gov- 
erning relations used to describe the viscous shock layer for a perfect gas. 

Boundary Conditions 

Conditions at the body surface .- The surface boundary conditions in terms of trans- 
formed variables are 

u = v = 0 (A14) 

T w = H w = Constant (A15) 

Conditions at the shock .- The shock conditions are determined by solving equa- 
tions (15) to (19). The transformed shock conditions become 

u = T = H= v = P = p = p=a.Q = a^ = a2 = l (A16) 


at T) = 1. 

Stagnation-Streamline Equations 

When downstream numerical solutions are required, it is necessary to have an 
accurate solution for the flow along the stagnation streamline. A truncated series, which 
has the same form as that used by Kao in reference 17, is used to develop the stagnation- 
streamline equations. The flow is assumed to be laminar along the stagnation streamline. 
The flow variables are expanded about the axis of symmetry with respect to the nondimen- 
sional distance | along the body as follows: 

P (£,??) = Pjfa) + P 2 ( t 7)£ 2 + • • • (A17a) 
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APPENDIX - Continued 


u(^,r?) = u 1 (t?)^ + . . . 

(A 17b) 

v(l,r?) = V!!??) + . . . 

(A 17 c) 

P(4 , 77 ) = P^t?) + . . • 

(A17d) 

TU, V ) = T^r?) + . . . - 

(A17e) 

h U,V) = h 1 (??) + . . . 

(A17f) 

= p 1 (7?) + . . . 

(A17g) 


The shock standoff distance is written as 


-S = n l,s + n 2,s« + - • ■ 


(A 18) 


Furthermore, £ is small and the curvature k is approximately 1 in the stagnation 
region. Consequently, the geometric relations (see fig. 1), including terms of order £, 
can be written as 


/3 ~ £ 


(A19) 


and 


a 



2n 


2 


il + n 


l,s 


Therefore, 

sin (a + (3) ~ 1 


and 


cos {a + 0) ~ 


2n 


2 i s_ 


1 + n 


l,s 


(A20) 


(A21) 


(A22) 
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The shock relations (eqs. (15) to (19)) in terms of expanded variables become 


M J(y - 1) + 2 


v s = v l,s + ' ' ' 


P l,s M oc 2 (y + 1) 


(A23) 


u s = u l,s^ + • • 


*4 


2n 


1 - 


1 + n 


1-_L 


l,s\ H 1,S ) 


(A24) 


Ps = Pl fS + p 2,s^ + 


+ 1 


yM 0 


l,s/ 


1- 


2n, 


2,s 


1 + n, 


M/\ i+ “l,s 


(A25) 


and 


T s = T l,s + ' • 




! ( r ~ ' ‘) 


4i- 1 


2 P 


l,s/ 


(A26) 


An examination of these equations shows that the equations for u g and p g con- 
tain n .2 g . This term cannot be determined from the stagnation solutions, since it is a 
function of the downstream flow. Consequently, a value must be assumed for n 2 g . In 
this study, it is assumed to be zero to start the solution, but this assumption is then 
removed by iterating on the solution with the previous shock standoff distances used to 
define n 2 jS . The effect of the downstream shock shape on the stagnation-point solution 
is elliptic rather than parabolic. 

Along the stagnation streamline, the second-order differential equations are written 
as 

^-+ a-, — + a 9 W + a, = 0 (A27) 

*, 2 1 2 3 

The coefficients are defined as 


s- momentum, W = u: 


a 


1 


1 ±1 1 (i , i) p i.s -iwi.. y 1 

5 1 <”> 1+n l,s" A* ljS 


(A28a) 
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a, 


_ n l,s 

1 dp l l 

n l,s 

1 + n l s n 

Mi d? 1 

1 + n 1;g r? 




o^fl 


l,s 


Mi 


°Vs 


"i 


- 2 Pl,s n l,s 

/V — * * 

4 t v>. 

n 2)S v dp 2 \ 

p l,s( 1+n l,s p ) u l,s p 

\ 2 P M . 

n i, s d7 ?y 

Energy (enthalpy), W = H: 




(A28b) 

(A28c) 


dp . 


dN 


a l = - 


Pr 


,1 ^ + ^ n l,s _ n l,s p l,s( N Pr,l) s v l,s p i N p r ,1 V 1 '£ 


Mi dr? N prjl 


d r\ 


1 + n 


l.s 


a2p l,s 


(A29a) 


P 1 



(A29b) 



(A29c) 


The preceding energy equation is for the thin-shock-layer approximation. When equa- 
tion (7) is used for the n-momentum equation, the following term must be added to equa- 
tion (A29c): 



N Pr,l^l dp 
Mj d V 


The remaining equations are written as follows: 
Global continuity: 


d_ 

dr] 


0 + n l,s r ?) J+lp l,s v l,s p in = - ( j + J ) n l ,s i 1 + n l,s^) 3p l,s u l,s p i a i (A3 °) 
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n- momentum: 


dp x 

dr] 


v l,s p l.s 


PiV 


Pi. 


ri 


dv-^ 

d?7 


(A3 la) 


When the thin-shock-layer approximation is made, the n-momentum equation becomes 


dp. 

— - = 0 
drj 


The Pg term that appears in equation (A28c) can be expressed as 


(A3 lb) 


dp 2 _ p l,s u l,s n l,s p l u l“ . 2p l,s u l,s n 2,s v l,s p l u l^ dv l . P2,s p l,s v l,s - dv l 

■ — - * ■' n * ■ - 1 ~ ■■ + . ■ 2 — . i - - — + ' / — / , U . V ■» 


- - 2 


2 P 


dv-i 


dv. 


drj 


p l 


,s 


1 + rjn x 


,s 


p l. 


1 + rjnj g dr] 


p l,s 


p l v l 


dr] 

(A31c) 


For the thin-shock-layer approximation, this term is 


#2 

dr] 



p i a i 2 
1 + 7 ? n l,s 


(A3 Id) 


These equations along with the equation of state constitute the nonlinear ordinary differ- 
ential equations that are solved along the stagnation streamline. 
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Figure 3.- Flow chart for solution sequence of viscous-shock-layer equations. 
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Comparison of eddy-viscosity profiles. Figure 9.- Comparison of tangential velocity profiles at 
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